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In this paper we present a generalization of the Goryachev - Chaplygin integrable case on a bundle of Poisson brackets, 
and on Sokolov terms in his new integrable case of Kirchhoff equations. We also present a new analogous integrable case 
for the quaternion form of rigid body dynamics' equations. This form of equations is recently developed and we can use 
it for the description of rigid body motions in specific force fields, and for the study of different problems of quantum 
mechanics. In addition we present new invariant relations in the considered problems. 

> 

. In this paper we present a generalization of the Goryachev -Chaplygin integrable case on a bundle 

! of Poisson brackets including (co)algebras so(4) and so(3, 1), and on the quaternion form of Euler- 
Poisson equations. Note that the generalization on the bundle is connected with the introduction of 
■ variables on the bundle analogous to Andoyer - Deprit variables. These variables are separating for all 
members of the bundle. We also obtain the L — ^4-pair of the generalized Goryachev -Chaplygin case 
of quaternion equations. In a particular case this pair is reduced to the L — A pair of the classical 
Ch ■ Goryachev -Chaplygin case. 



Let us consider a generalization of the classical Goryachev -Chaplygin case of Euler- Poisson 



^ ' equations 

U\ M = Mx|i+ 7 xM j = Jx dH 



dM 1 cV ' ' dM 

for the zero value of the area integral (M, 7) = 0. We add a gyrostatic moment and singular term to 
the Hamiltonian 

H = \{Ml + Ml + 4M|) + AM 3 + /x 7 i + \\, (1) 
where A, /x, a = const. The additional integral is of order three with respect to the moments 



M 3 + |) (Ml + Ml + -|) - M M l73 . 

D. N. Goryachev himself in the paper [2] showed the generalization on zero value surface of gyro- 
static moment A = (for A 7^ 0, a = it was shown by L. N. Sretensky 9 ). The complete form of 
generalization Q was considered by I. V. Komarov and V. B. Kuznetsov jSj. They also presented some 
quantum mechanical interpretation of the singular term. 
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The Poisson bracket of variables M,-f is defined by algebra e(3) and has the following form 

{M i ,M j } = -e ljk M k , {M i)7j } = -e iifc Tfc, H,^} = 0. (2) 

Let us present the generalization of the Goryachev-Chaplygin case on the bundle of brackets 
of the form 

{Mi, Mj} = -£ij k M k , {Mi,7j} = -e ijk >y k , {ji,^} = -xe ijk M k . (3) 

At x = 1 these commutation relations correspond to the algebra so(4). 
We can present the Casimir functions of bracket in the form 

Fi = (M, 7 ), F 2 = 7 2 + xM 2 . (4) 

Now let us construct the variables on bracket analogous to the Andoyer-Deprit variables pQ. 
Suppose the component M3 is equal to the momentum 

L = M 3 . (5) 

The variable / canonically conjugate to the variable L ({I, L} = 1) on subalgebra so(3) with the 
generators M\ , M2 , M3 is constructed by integration of the Hamiltonian flow with Hamiltonian func- 
tion J$? = L 

^ = { Ml ,L} = M 2 , ^ = {M 2 , L} = -Mi, 

(6) 

if = (M 3 , L, = 0. 
Hence using the commutation relation {Mi, M 2 } = — M3 we obtain 



Mi = VC 2 -L 2 sinZ, M 2 = \J G 2 - L 2 cos/, (7) 

where G 2 = M 2 + Af| + M 2 is the Casimir function of subalgebra so(3). 

Suppose G is the second momentum and construct the canonically conjugate variable g. We 
chose H = G as a new Hamiltonian, and the corresponding flow on the whole bundle Jz? x has the form 

dM _ n dj _ 1 

"AT ' d# ~~ G (8) 

where 5 is the variable canonically conjugate to G. 

According to 0, M does not depend on g, and using equations and Casimir functions 
we obtain for 7 the equations 

7 = -^-M + ^(M x e 3 sin g + GM x (M x e 3 )cossi), 
G 2 G 

c 2 -xG 2 -^ W 
G 2 -L 2 



« 2 = ^ 7T^> e 3 = (0, 0, 1), 



where c 2 = xiW 2 + 7 2 , and H is the traditional notation of the area integral H = (M, 7) = c±. 

Thus ©10)0 define the symplectic variables on the whole bundle ££ x corresponding at x = 0, 
c = 1 to the well-known Andoyer-Deprit variables in rigid body dynamics. 
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Using ((SJ), @ we obtain the generalization of the particular Goryachev-Chaplygin integrable 
case on the bundle J£ x . We chose the Hamiltonian in the form 

= I(G 2 + 3L 2 ) + AL + a(cos/cosc/ + ^sin/sing), (10) 

where a, A are constants. 

In (|1U|) we add the linear with respect to L term. It is interpreted on algebra e(3) as the component 
of the gyrostatic moment 

We can separate variables in system ()10|) . Indeed, let us apply the canonical change of variables 

L = pi+p 2 , G = p 1 -p 2 , qi = l+g, q 2 = l-g. (11) 
Now Hamiltonian (jlUj) can be presented in the form 

3 3 2 2 

= A (pi singi +p 2 sing 2 ). (12) 

Using (J5J), ((7J), © we present Hamiltonian (jlUj) as a function of variables iVf, 7 on the zero surface of 
area integral (M, 7) = H = 

^ = I(M 2 + M 2 2 + 4M 3 2 ) + AM 3 + /i-^. (13) 

2 It I 

The additional integral in this case has the form 

F = (M 3 + |) (M x 2 + Mf) - ixM x |^|. (14) 

On algebra e(3) we have I7I = 1 and obtain the classical Goryachev-Chaplygin integrable case. 
In the classical case this method of variables' separation was suggested by V. B. Kozlov [2]. 
Singular term (^Q) is generalized on the bundle in the following way: 



2 

H = \{Ml + Mf + 4M|) + AM 3 + fj,^L + h—Y, 
2 | 7 | 2 7 | 

F = (m 3 + |) (M 2 + Mf + aj|) - ^M l73 , 

although variables' change in this case doesn't produce the separation of variables (at least we 
don't know such separation). 

Recently Sokolov and Tsyganov present a new generalization of particular integrable case Q 
on bracket (J2J) with the Hamiltonian containing the quadratic cross terms with respect to M,~f. 

The most general form of the integrable family in this case is presented as the following Hamil- 
tonian 

H = i (m\ + Mf + 4Mf + 4) + AM 3 + nxl\ + M272+ 

2V 7; V (15) 

+ai(2M 3 7i - M173) + a 2 (2M 3 72 - M 273 ). 

The additional integral is 

F = (m 3 + 0171 + a 272 + |) (m? + Mf + -|) - (m M i + M2M 2 ) 73 . (16) 
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At e = we can separate variables in Hamiltonian (|15[) using the linear combination of Andoyer- 
Deprit variables analogous to Indeed, we can show that 

H = p - p2 (/(Pl>gl) -0CP2,«2)), 

/ = 2aipf + 2 (f + ai sin ^ + a^cosq^Pi + (Mi s ' m H + M2Cosgi)pi, (17) 

5 = 2aiP2 + 2 (| - a i sin 92 - a 2 cos g 2 ) p| ~ G"i sin 92 + A*2 cos <? 2 )P2 • 

As we note above, we don't know such separation at e 7^ 0. 

If the gravity field is absent (/xi = fi 2 = 0), then the integral is presented as a product of factors 
F = kik 2 , where 

h =a M 3 + a l7l + 0272 + ^7, fc 2 = M 1 2 + M| + ^. (18) 

2 73 

We can easily show that the equations of motion for k\ , k 2 have the form 

fci = 2(ai72 - 027i)fci> ^2 = -2(oi7 2 - a 2 -yi)k 2 . 

{M,y)=0 

Thus the equations /ci = and k 2 = define the invariant relations of system (|15|) at = fi 2 = (and 
fc 2 = is the particular invariant relation on zero level of (M,~f) = 0, and k\ = is general invariant 
relation). It seems that, invariant relations Q18JI and the additional integral in the form F = k\k 2 
where introduced by authors. 

Using the observation from the first section we can generalize integrable case ()15[) and invariant 
relations on bundle of brackets (jSJ). We shall make the following substitution in Hamiltonian (|15|). 
integral (jlfij) . and invariant relations (|18|) 



7i^rr, |7| = \/7i +72 + 7 3 - 
171 v 

In this case at e = we can also separate variables using the analog of Andoyer-Deprit variables and 
substitution 

We should note that unlike the generalization of Kovalevskaya case pQ, where the Hamiltonian 
and integral explicitly depend on the parameter of bundle x, in the generalization of Goryachev- 
Chaplygin case the integrals don't depend on the parameter of bundle. 

Let us consider another possibly not so natural case of motion equations of rigid body with 
a potential linear with respect to Rodrig- Hamilton parameters and not with respect to direction 
cosines P 

3 

H= i(AM,M)+^r i A il n = const, (19) 

i=0 

The equations of motion in this case are 

M - M v dH + 1 \ * dH + 1 dH \ 1 \ dH 
M " MX 9M + 2 AX 75A + 2^ A -2 Ao 75A' 

(2U) 



A °-"2v A '5M> A " 2 AX dM + 2 Xo dM- 



where A = (Ai, A2, A3 



24 



REGULAR AND CHAOTIC DYNAMICS, V. 7, J\U, 2002 



GENERALIZATION OF THE GORYACHEV CHAPLYGIN CASE 




These equations are Hamiltonian one with the Poisson bracket 

{M i ,M J } = -e ijk M k , {M, , A } = \ \ , 

(21) 

{Mi,Xi} = -^(sijk^k + <%A ), i,j,k = 1,2,3, 

corresponding to the Lie algebra so(3) (S> s M 4 . The relations between quaternions and Euler angles 
have the form 

. e i> + <p . . e ^-<p 

A = cos-cos — - — , Ai=sm-cos — ^ — > 

. . e . 4>-(f . e . jj + v 

A 2 =sm-sm — - — , A3 = cos ^ sin — ^ — • 

In the classical mechanics such potential aren't presented because its dependance on the body position 
is nonunique (more exactly it is double- valued function). To justify the study of such equations we can 
refer to the problems of quantum mechanics and point masses dynamics in curved space S" 3 , and on 
some formal technics of L — A-pair construction . It turns out that after the reduction of order 
of system (|19[) we obtain the general Euler -Poisson equations with additional terms having different 
physical interpretations 

As an interesting feature of system ()19|) we note that using linear with respect to Aj transforma- 
tions we can reduce the general form of the potential 

3 

V = Y,n\i (22) 

to the form 

V = r A . (23) 

Indeed, the linear transformations of quaternion space Aj (preserving commutation relations and 
quaternion norm) of the form 

Ao = i? _1 (r A + riAi + r 2 A 2 + r 3 A 3 ), 
Ai = i? _1 (r Ai - nX - r 2 A 3 + r 3 A 2 ), 

A 2 = i?~ 1 (r A 2 + riA 3 - r 2 A - r 3 Ai), (24) 
A 3 = i? _1 (r A 3 - r x A 2 + r 2 Ai - r 3 A ), 
R 2 = rl + r\ + r\ + r\ 

reduce potential ()22fi to form (|23|) . The existence of such linear transformation is the remarkable 
property of quaternion variables and bracket (|21j) . the analogous transformation does not exists for 
the brackets of algebras e(3) and so(4). 

In general dynamically unsymmetrical case a\ ^ a 2 7^ a 3 ^ a\ system (|19|) seems to be noninte- 
grable and none of two necessary additional integrals exists. However, this was not proved, and the 
proof by various reasons is not simple. Note that even the application of the Kovalevskaya method for 
system (|19j) is not quite analogous to the classical Euler -Poisson problem. Generally speaking even in 
the Euler - Poinsot case the solution branches on the complex plane of time (with the exponent 1/2). 
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At a\ = a 2 the linear integral always exists 



Fx = M z {rl + r\ + rl + rf) + N 3 (rj + r\ - r\ - rf )+ 
+ 2N 2 (rir - r 3 r 2 ) - 2Aq(?qr 2 - r r 3 ), 



(25) 



where iVj are projections of kinetic moment on the fixed axes. For potential Q23|) this integral has the 
natural form 



It turns out that (linear) integral ()26j) corresponds to the cyclic variable ip + ip P . Rauss reduction 
with respect to this cyclic variable results in Hamiltonian system on algebra e(3) on zero surface of 
area integral (M,~f) = with the Hamiltonian 



where c is a constant value of integral (|26|) . Hamiltonian (|27|) corresponds to the addition of a 
gyrostatic term linear with respect to M and singular term into the general Euler - Poisson equations. 
The physical meaning of singular term is discussed in 

Integrable cases of system (|19|) (equivalent to integrable cases of system (|27|)) are presented in PQ. 
Here we just present the generalization of Goryachev- Chap ly gin case. 

Hamiltonian and additional integral are 



Integral F% commute with integral F\. Under the reduction to system (|27|) this case becomes a member 
of generalized family (fT|). 

Remark 1 . Adding a constant gyrostatic moment along the axis of dynamical symmetry in l|28[) we obtain 
an integrable case corresponding to the generalized Sretensky case in Euler -Poisson equations. Additional 
integrals can be easily obtained with the help of the lifting procedure described in . 

Note that the indicated "Goryachev -Chaplygin case" for the quaternion Euler - Poisson equa- 
tions is the general integrable easel Thus we can use it for some algebraic constructions (L — A-pair 
constructions and all that) described below. 

Lax representation of systems on algebra su(2, 1). We start from some formal method 
of L-A pair construction based on existence of two compatible Poisson brackets for the same system 
of Hamiltonian equations on Lie algebra; in this case the system is called bihamiltonian. The detailed 
exposition of this method and its applications to various problems of rigid body dynamics are presented 



F x = M 3 - AT 3 . 



(26) 



H = \{Ml + Mf + a 3 M|) + c(a 3 - 1)M 3 + r o72 + \^, 
2 2 7 | 



(27) 



H = ±(M? + Mf + 4M|) + r A , 

F 2 = M 3 (Mf + Mf) + r (M 2 Ai - MiA 2 ). 



(28) 



in 0E]. 
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Let's consider space of complex matrixes 3x3 with the basis 



Mi = 


f° ¥ 

¥ 





M 2 = 


( ° 

1 

2 










^ o 



0/ 





(-¥ 


\ 









M 3 = 


¥ 












0/ 



M 4 



/-ii o 
6 

-¥ 


\ 




\ 


¥) 





/ 


i> 















2 










Pi = 







P 2 = 











{-¥ 


°y 




[-Ixi 





) 




/ 


°\ 








o \ 









1 













P 3 = 




2 


P 4 = 






-¥ 






[ixi 


°/ 






J 





(29) 



With respect to the standard matrix commutator [■ •] such matrices generate a semi-simple algebra 
with Cartan decomposition Jz? = H + V, where subalgebra H = su{2) © su(l) is generated by 
matrices Mj, and V = C 2 is generated by matrices Pj. Here x is a parameter determining some 
bundle of algebras linearly dependent on x. For x > these algebras are isomorphic to algebra su(3), 
for x < they are isomorphic to algebra su(2, 1), at x = they are isomorphic to the semidirect sum 
(so(2) © su(l)) ® s C 2 . 

Using the semi-simplicity we can identify algebra with coalgebra by means of the inner product 
(Killingform) 

5 = -Tr(X-Y), X, Y€ if. (30) 

Let's denote coordinates in coalgebra as mi, m2, m.3, m4, pi, P2, P3 ; P4, then after the identifica- 
tion we obtain a matrix (an element of algebra): 



X 



/-i(m 3 + m 4 ) imi+m 2 ^(pi 
imi - m 2 i(m 3 - m 4 ) -(p 3 



IP2) 



m) 



V 



-Pi - W2 



~P3 ~ WA 



The corresponding Lie-Poisson bracket, more precisely, the bundle of brackets linearly dependent 
on parameter x has the following form for the coordinate functions of coalgebra 



{mi, rrij} = e ijk m k , {mi, m 4 } 

{rrn, pj} = ^(e ijk pk - SijPi), {m, Pa} 

{Pi, Pj} = \x{e ijk m k + 3£y 3 m 4 ), {pu Pa} 



0, 
1 



1 



j Ji k — 1, 2, 3 

x(rrii - 3(5 i3 m 4 ). 



(31) 



In correspondence with the general method developed in |21E| , we shall construct L matrix, and its 
invariants define a commutative set for the whole family of brackets { }g + A({ }\ + { } a ), where a £V 
is a shift of the argument. The bracket {}# is different from (|31|) ; in this bracket variables pi are 
pairwise commuting. Let's assume x = 1 restricting ourselves to the problem with a real dynamical 
sense. Then we can present L-matrix as: 



L = (hA + v + aA 2 ), 
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where 





^-7(7773+7774) 


imi+m,2 




H 


im\—rrL2 


7(7773—7774) 


w 




\ 


277774/ 



( ' 


Pi 

P3 


- «P2\ 

- «P4 


\-Pl 




~P3 






o / 











ai 


- m 2 \ 


(~ai 






a 3 


— ia4 ] 




-ia 2 


-a 3 


— ia 4 


/ 



(32) 



Among invariants of matrix L under arbitrary shift there is a linear on m^ function of the form 

F 1 = m 4 a 2 + (m, -y a ), (33) 



where a 2 = a 2 , m = (mi, m 2 , 013), and the components of vector 7 a are the functions of a, of the 

i=l 

form 

7 a = (2(oi0 3 + a 2 a 4 ), 2(-a 2 a 3 + aia^), af + af - a\ - a 2 ). 
Let's suppose a\ = a 2 = 0. In this case integral l|3*37) has the form 

F\ = m 3 + m 4 . (34) 



The following square-law invariant of matrix L, we will choose as a Hamiltonian 

O \ OOO o 

F 2 = Tr(h + 2va) = m 1 + m 2 + m 3 + 3m 4 + 04/34 + 03/73. 



(35) 



It defines some (formal) integrable Hamiltonian system on the family of brackets {}q + A({}a + {}a)- 
We can easily find matrix A for this system [5]. 

The reduced system and nonlinear Poisson structure. In order to proceed from the 
found formal system to the generalization of Goryachev-Chaplygin case [7j, we shall make the re- 
duction using linear integral (|34j) directly in obtained L and A matrices. In the general case this 
reduction is not Poisson reduction [5]. We shall make the following substitution in matrix L (|3Uj) and 
in Hamiltonian 1)35(1 

777,4 = — 777,3 + c, c = const. 

We obtain L matrix and Hamiltonian of integrable system on subalgebra mi, m 2 , 7713, pi, /7 2 , 733, p± 
with gyrostat, the gyrostatic moment being equal to c: 



/ — i\c (imi+m 2 )X Pi~ip2 \ 

(imi— m 2 )A i\(2ni3—c) P3—ip4+(a3+ia,4)\ 2 

\ -pi-ip2 -P3-^4-(fl3+04)A 2 -2i(m 3 -c)X / 



H = m 2 + m 2 + 4m 3 — 6771,3c + 204/74 + 203/73 



(36) 



(37) 



^(47773 — 3c) —77771 — m 2 

A = dH = I —imi + m 2 —7(47773 ~~ 3c) 

03 + 704 
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Let's consider Poisson structure, defined by a bracket {}#. The corresponding Hamiltonian system 
with Hamiltonian Q37|) has another linear integral 



^3 = ^-3- (m, 7) 
7 = (2(f2P4 +P1P3), 2{P2P3-Pm), P3+PI-P1-P2)), 



, 2 2 2nn ( 38 ) 



using this integral we can carry out a usual procedure of order reduction (such as Routh reduction or 
reduction on moment). We can carry out the reduction in the simpler way using algebraic form, if we 
choose as new variables the integrals of vector field v = {•, F3} of the form [3] 

M lPl + M 2 p 2 M 2 pi - M lP2 

Ai = - — , Ao = - — , A3 = Ivli, 

VP1+P2 VP1+P2 (39) 



si=Ps, s 2 =P4, s 3 = ±ypl +p 2 , 

having the nonlinear commutation. We describe the transformation (|39[) in |Hj 

p 

{Ki, Kj} = EijkKk + £ij3 — 

s 3 



4 (40) 



where F3 = (X, s)«3 is a Casimir function of structure (|4()|l. 

At zero value of "the area integral" , i. e. for F3 = s)s3 = bracket ()4()|) coincides with a usual 
Lie -Poisson bracket, defined by algebra e(3) = so(3) © s M 3 , and Hamiltonian (|3*7|) in variables P9*|) 
coincides with the Hamiltonian of the classical Goryachev-Chaplygin case [Zj 

H* = \H = ±(Kl + K% + 4Ki) - a 4 s 2 - a 3Sl , 03, a 4 = const. (41) 

Goryachev case with a singular term. We can remove the nonlinearity of structure (|4U|). 
appearing for F3 7^ on the fixed level .F3 = c using the transformation 

L = K-cf 3 . 

After the transformation bracket (|4()j) gets the form of algebra e(3), and Hamiltonian H* gets the form 

H* = \(L\ + L 2 2 + 4L§) - a A s 2 - a 3 si + 3cL 3 + \^ (42) 

Hamiltonian (|42j) can be interpreted as some generalization of the Goryachev -Chaplygin case, at 
(L, s) = 0, with the additional terms, linear on L3 and corresponding to the gyrostatic moment. The 
integrable generalization with the gyrostatic moment only was described by L. N. Sretensky in [H], the 
generalization with singular potential only was presented by D. N. Goryachev himself 0, the general 
case, when it is possible to add to the Hamiltonian both terms with arbitrary independent coefficients, 
was described in paper [Sj. 

Thus, the presented L-A pair is valid for the generalizations of the Goryachev -Chaplygin case. 
It is different from the mysterious L-A pair described in paper [3] which is obtained by deletion of a 
row and column from the relevant pair of the Kovalevskaya case. 

The authors thank A. V. Bolsinov, A. V. Tsyganov, V. V. Sokolov for the valuable remarks and 
numerous discussions. 
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